Conductance as a Function of the Temperature in the Double Exchange Model 
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We have used the Kubo formula to calculate the temper- 
ature dependence of the electrical conductance of the double 
exchange Hamiltonian. We average the conductance over an 
statistical ensemble of clusters, which are obtained by per- 
forming Monte Carlo simulations on the classical spin orien- 
tation of the double exchange Hamiltonian. We find that for 
electron concentrations bigger than 0.1, the system is metal- 
hc at all temperatures. In particular it is not observed any 
change in the temperature dependence of the resistivity near 
the magnetical critical temperature. The calculated resistiv- 
ity near Tc is around ten times smaller than the experimental 
value. We conclude that the double exchange model is not 
able to explain the metal to insulator transition which exper- 
imentally occurs at temperatures near the magnetic critical 
temperature. 

PACS number 71.10.-w, 75.10.-b 



I. INTRODUCTION 

Materials that present extremely large magnetoresis- 
tance have potential technolog- 

ical applications. As mixed valence compounds of the 
form Lalt^Al+ Mnlt^Mn^+Ol-- (where A can be Ca, 
Sr or Ba) show colossal magnetoresistance there is a re- 
newed interest in these oxides with perovskite structure 
There is a big correlation between the magnetic 
and the transport properties of these oxide manganites. 
For 0.1 < X < 0.5 and low temperatures, the system is 
metallic and presents ferromagnetic order. As the tem- 
perature (T) increases the system becomes insulator and 
paramagnetic. The metallic (insulator) behavior is de- 
fined in the sense that dp/dT > {dp/dT < 0), being p 
the electrical resistivity. In the insulating phase the resis- 
tivity is bigger than the Mott resistivity ^ lOOO/ifi — cm 
1^. The magnetic transition occurs at a x-dependent 
critical temperature Tc ~ 300K. The metal to insulator 
transition occurs at a temperature very close to Tc. For 
X and low temperatures the system is a layer antifer- 
romagnet with ferromagnetic coupling inside planes. At 
X < 0.1 phase separation between hole-rich and hole-poor 
regions has been predicted The electron-electron 

interaction plays an important role at low temperatures, 
and it produces charge ordering at particular values of 
X < 0.5 and always for values of x > 0.5 0|. 

In the Mn oxides the electronically active orbitals are 
believed to be the Mn d orbitals, and the mean d occu- 



pancy is A — x. The Mn ions are located at the corners of 
a simple cubic lattice, and they feel the cubic crystal sym- 
metry, which splits the d orbitals into a t2g triplet and a 
eg doublet. There is also a strong ferromagnetic Hund's 
rule coupling which align all electron spins in the Mn 
d orbitals. The physical picture is that three electrons 
fill up the t2g levels forming a core spin S of magnitude 
3/2 and the rest of the electrons go to the Cg orbitals. 
For small values of x, the perovskites show a long-range 
Jahn- Teller order which selects a preferred combination 
of the Cg orbitals and therefore it is possible to assume 
that the electrons move only through one d orbital. 

To explain the ferromagnetism in these materials Zener 
1^ introduced a double exchange (DE) mechanism, in 
which the electrons get mobility between the Mn ions 
using the magnetically inert oxygen as an intermediate. 
The oxygens are located at the center of the lines connect- 
ing the Mn ions. This conduction process is proportional 
to the electron transfer integral and due to the strong fer- 
romagnetic Hund's rule coupling it is maximum when the 
two cores spins involved in the process are parallel and 
it is zero when they are antiparallel. So in the DE model 
ferromagnetic coupling between Mn^^ and Mn'^'^ arises 
from the hopping of the electrons in the Cg orbitals. Be- 
cause the alignment of spins on neighboring sites favors 
electronic motion, the ferromagnetic ground state maxi- 
mize the electron kinetic energy. When the temperature 
increases, the DE model undergoes a phase transition 
towards a paramagnetic state. In this phase the core 
spins are randomly oriented and fluctuate at frequencies 
related only to the temperature. In the paramagnetic 
phase the electrons minimize their kinetic energy. The 
DE model was more precisely formulated and extended 
by Anderson and Hasegawa |^ and by de Gennes [ p^ . 

Mean field theories have shown that the ferro to para- 
magnetic phase transitions at T,. is accompanied by a 
change in the temperature dependence of the resistivity 
|pd]-p^. However these calculations do not show metal 
to insulator transition at temperatures near T^. These 
mean field calculations do not take into account the ef- 
fect of the Berry phase arising from particle motion in 
a spin background |l^ . Also the possible localization of 
electrons by the spin disorder in the paramagnetic phase 
is neglected. Varma proposed that random hopping 
in the paramagnetic phase is sufficient to localize elec- 
trons and induce a metal to insulator transition at the 
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magnetic critical temperature. However Li et al ||17|] have 
studied the mobihty edge of the DE model in the limit 
of r — > oo and they have found that random hopping 
alone is not enough to induce Anderson localization at 
the Fermi level in the range 0.1 < x < 0.4. Similar re- 
sults were reported by AUub and Alascio ||T^ using the 
Ziman criterion [ p9[ . 

In this work we study the temperature dependence of 
the electric conductance of the DE model for different 
values of the electron concentration. We find that the 
system s metallic at all temperatures and, contrary to 
the mean field results, we do not observe any feature in 
the temperature dependence of the resistivity near the 
magnetic critical temperature. 

We calculate the dc conductance, G, by using the Kubo 
formula. We perform Monte Carlo simulations on the 
classical spin orientations of the DE model |2^] , and for 
each electron concentration and temperature we obtain 
a statistical ensemble of clusters of Mn ions. We obtain 
the conductance by averaging over several configurations. 
The temperatures we are interested are T < 500A' (Tc ^ 
300K) and these temperatures are much smaller than 
the electron Fermi temperature for x > 0.1. Therefore 
we always consider, both in the Monte Carlo simulations 
and in the conductance calculations, that the electron 
temperature is zero. 

In our model the Mn ion spins are treated as classi- 
cal and quantum effects, as absorption and emission of 
spin density waves, are not correctly described. Although 
quantum effects should be important at very low temper- 
atures | pl| , they do not matter at temperatures near Tc- 

The paper is organized as follows, in section II we intro- 
duce the DE Hamiltonian and describe the Monte Carlo 
calculations. In section III we give some remarks about 
the use of the Kubo formula in the present context. Sec- 
tion IV is dedicated to present and discuss the results of 
the paper, and we finish in section V with a summary. 

II. DOUBLE EXCHANGE HAMILTONIAN. 

For large Jahn- Teller splittings, the electronic and 
magnetic properties of the Mn oxides are described by 
the following ferromagnetic Kondo lattice Hamiltonian 
§, 

<i,j>,a 

^ Jh ^ C^^a '^a-(y' ^i.c' ■ I (1) 

where C^^ creates an electron at site i and spin cr, Sj 
represents the classical core spin at site i, t is the hopping 
amplitude between nearest-neighbor sites and Jh is the 
Hund's rule coupling energy. In the limit of infinite Jh, 
Eq.(l) becomes the DE Hamiltonian, 



Hde = -Y1 {u,,C+C, + h.c}j , (2) 

Here creates an electron at site i with spin parallel 
to Si, and the hopping amplitude acquires a Berry phase 
and it becomes a complex number given by | [l5|] , 

f,^.=i(^cos|cos|-Hsin|sin|e^W'-^^)) , (3) 

where 9i and are the angles which characterize the 
orientation of S^. This complex hopping appears after 
rotating the conduction electron spins so that the spin 
quantization axis at site i is parallel to Si, and then 
project onto the spin parallel to Sj. 

In order to calculate the electrical conductance, for 
each T and x, it is necessary a statistical ensemble of 
clusters, which simulates the thermal fluctuations. Each 
cluster is characterized by a set of core spin {Si} and its 
chemical potential ii. We obtain the clusters by perform- 
ing Monte Carlo simulations on the variables 6i and (jn 
of the Hamiltonian Eq.(2). In each Monte Carlo step it 
is necessary to diagonalize a matrix of size equal to the 
number of Mn ions in the clusters and the DE energy 
is the sum of the eigenenergies of the occupied electron 
levels. In this process we assume that the electron Fermi 
energy is much bigger than the temperatures of interest in 
this work (at x ~0.2, Tperrm ^ 2000ii: and ~ 300K). 
The diagonalization imposes a restriction on the dimen- 
sion of the unit cell used in the simulations. Recently two 
of us have studied, by using Monte Carlo simulations, the 
magnetic phase diagram of the DE model . To avoid 
finite size problems in the simulation, the following ap- 
proximation for the energy has been obtained and used 
for the DE energy, 

E ^ -2t{d+d,)o E COS ^ - a2^ii-U,,f . (4) 

(id) («J'> 

Here (C+Cj)o and 02 are quantities which depend on x 
and do not depend on temperature, t is the average of 
the absolute value of the hopping amplitude and 9ij rep- 
resents the angle formed by the core spins located at sites 
i and j. In reference [20[ it was showed that expression 
(4) is a good approximation to the DE kinetic energy. In 
that work it was also concluded that the magnetic criti- 
cal temperatures of the DE model are in the range of the 
experimental ones, and it was obtained that the complex 
phase of the hopping amplitude has a negligible effect in 
the value of Tc. 

Using the conclusions of the cited work, in this paper 
we perform Monte Carlo simulations on the variables 9i 
and (/)i, using the expression (4) for the DE energy. The 
simulations are performed in N x N x N cubic lattices. 
Technical details about the Monte Carlo calculations are 
given in reference |2^, here just say that typically 5000- 
7000 Monte Carlo steps per spin are used for thermal 
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equilibration. The different clusters which form the sta- 
tistical ensemble are chosen every 100 steps per spin after 
equilibration. 



III. D.C. CONDUCTANCE VIA KUBO FORMULA 

In calculating the conductance of the system we use 
the standard Kubo formula |^,|2|. The static electrical 
conductivity at chemical potential fi is given by 



;I - H - Si(fe) - S,(£)]G(e) = I 



(7) 



G = a,,(0) = -2-Tr {hv,)lmg{ti)ihv,)lmg{fi) 
h L 



(5) 



where Im Q (/i) is calculated from the advanced and re- 
tarded Green functions 

Img(M) = ^p^(/.)-^-^(Ai) 

and the velocity (current) operator Vz is related to the 
position operator z through the equation of motion 



H,z 



(6) 



H being the Hamiltonian in Eq. (2). 

Numerical calculations are carried out for a bar ge- 
ometry connecting the N x N x N cluster to two ideal 
semiinfinite leads oi N x N cross section. This provides 
complex selfenergies at opposite sides of the sample . 
They are first calculated for the normal modes of the 
lead and then transformed to the local tight-binding ba- 
sis. Specifically, the retarded selfenergy due to the mode 
of wavevector (fcj,, fcj,) at energy e is given by: 



^{kx,ky) = ^ - e{kx,ky) - At'^ ~ (e - £{kx,ky))' 
within its band and by: 



Y.{kx,ky) = ^ ( e - e{kx,ky) T \/ (e - e{kx, ky)Y - 



outside the band (minus sign for e > £{kx, ky) and plus 
sign for e < £{kx,ky)), where e{kx,ky) = 2t{cos{kx) + 
cos{ky)) is the eigenenergy of the {kx,ky) mode). The 
transformation from normal modes to the local tight- 
binding basis is obtained from the amplitudes of the nor- 
mal modes: 

2 

< {nx,ny)\{kx,ky) ^ ^ ^ sin(fca;n^) sm{kyny) , 

where and Uy represent the tight-binding orbital posi- 
tion. Once the selfenergy matrices introduced by the left 
(right) semiinfinite leads are determined, the retarded 
Green function matrix of the sample G(£) is defined by 
the following set oi N x N x N linear equations: 



where Si(i.)(i?) stand for the selfenergy matrices intro- 
duced by the left (right) semiinfinite leads. This set of 
equations is efficiently solved using a layer by layer inver- 
sion scheme that takes advantage of the band structure 
of the coefficients matrix. The advanced Green function 
matrix is simply the conjugate of the transpose of the 
retarded one. 

The last ingredient that is necessary for the evalua- 
tion of the Kubo formula is the velocity operator. It is 
obtained through Eq.(^) once the position operator z is 
known. This operator is determined by the spatial shape 
of the electric potential energy. Taking advantage from 
the fact that the detailed form of the electric field does 
not matter within one-electron linear response theory, an 
abrupt potential drop at one of the two cluster sides pro- 
vides the simplest numerical implementation of the Kubo 
formula Certainly, Eq.(|^) shows that nonvanishing 
elements of the velocity operator are restricted to the 
two layers at the sides of the potential drop. Further- 
more, the trace appearing in Eq.(|^) makes the knowl- 
edge of the Green functions on the same restricted set 
of sites enough for the evaluation of the conductance. 
Consequently, Green functions are just evaluated for two 
consecutive layers at the cluster boundary. 



IV. RESULTS 

We start showing the results corresponding to a perfect 
ferromagnetic system. In Fig. la we plot the conductance 
as a function of the chemical potential for a cluster of 
size 7V=20 at T=0. The conductance is finite for chem- 
ical potentials in the range —6t < fJ, < Qt. This is the 
range of energies where the density of states of the perfect 
' system is finite, see Fig. lb. Since at T=0 all core spins 
are aligned, the transport in the system is ballistic and 
the conductance is just limited by the size of the cluster. 
Therefore, at T=0 the conductance is proportional to the 
, number of transport channels at the Fermi energy, which 
for large clusters increases as N'^. We have verified this 
behavior for TV > 10 and also we have checked that for 
system sizes bigger than N=10, the overall shape of G 
does not depend on the cluster size. 

In the opposite limit we have calculated the conduc- 
tance in the paramagnetic phase at T — > cx). In this limit 
the orientation of the local spins {S^} is random. We 
choose Si to be uniformly distributed on a sphere, i.e. 
the probability P of having a core spin with azimuthal 
angle 0^ is l/27r, whereas the polar angle distribution 
verifies P{cos9i)=l/2. 
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FIG. 1. a) Energy dependence of the conductance for the 
T = case (empty squares), the T — > oo case (filled dots), 
and the T — > oo case without Berry's phase (empty dots), b) 
Density f states as a function of the energy for the T = case 
(empty squares) and the T —f oo case (filled dots). 

In Fig.l we plot the conductance and the density of 
states as a function of the energy in the T — > oo limit. 
This quantities are averaged on many configurations of 
clusters of size N=20. The density of states shows an ef- 
fective band edge at Eb ~ — 4t. This occurs because the 
average value of the absolute value of the hopping ampli- 
tude at T ^ cx) is < \tij\ >=2/3i. Although the band 
edge starts at —At, G is different from zero only for values 
of the chemical potential in the range — 3.6t < ^jl < 3.6<. 
The difference between the band edge energy and the 
minimum energy with G ^ is due to the fact that all 
the states with energy —At < E < —3.6t are localized 
and do not contribute to G. Therefore Ec ~-3.6i is the 
T — > cx) mobility edge of the DE model. This result is 
in agreement with localization length calculations [ p7[ , 
which located the mobility edge at \Eb\ ~ 3.56i. Our 



calculations show that \Ei,\ is constant for values of N 
bigger than N=10. This implies that for energies lower 
than Eb the localization length is shorter than 10 lattice 
parameters. For N >10, the overall shape of the average 
conductance as a function of energy is almost indepen- 
dent of the cluster size, and its value increases linearly 
with TV. This implies that in the T ^ oo limit the sys- 
tem is metallic, the electron transport is diffusive and the 
system verifies the Ohm's law: the cluster resistance is 
proportional to the cluster length and inversally propor- 
tional to the cross-sectional area. The proportionality 
constant is the resistivity p. 
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FIG. 2. Conductance as a function of the cluster size, for 
X = 0.5 and T — > oo. The dashed line is only a guide to the 
eye. 

As an example of this behavior we plot in Fig. 2 the 
average conductance at /.t = as a function of the cluster 
size N . The linear behavior of G versus N down to A^=4, 
implies that the elastic mean free path is smaller than 4 
lattice parameters. The slope of the straight line is the 
conductivity, in lattice parameter units, of the DE model 
at T — > cxD and a;=0.5. 

In Fig.l we also plot the average conductance in the 
T ^ oo limit of the DE model but neglecting the Berry's 
phase in the hopping amplitude i.e. ti^ —^ \tij\. In this 
case G is bigger than in the case of complex hopping. 
Also it seems that the mobility edge appears at lower 
energies than in the complex hopping DE model. 

The previous calculations correspond to the T=0 and 
T — > (X) limits. Now we present the results for G as a 
function of T, for different values of a;. In these case the 
average of G is done with an statistical ensemble of clus- 
ters which are obtained from Monte Carlo calculations 
as described in Section II. When T increases the disorder 
in the system increases and two main effects occurs: a) a 
T-dependent mobility edge appears at low energies and 
b) the extended states acquire an elastic mean free path. 
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£, which decreases with T. 
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FIG. 3. Temperature dependence of the conductivity for 
the case of a; = 0.5 and A'^ = 16 and A'^ = 24. The magnetic 
critical temperature Tc is pointed with an arrow. Above this 
temperature conductivity is independent of A^. 
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FIG. 4. Temperature dependence of the conductivity for 
the case of a; = 0.2 and A^ = 16 and A^ = 24. 

With respect to point a) we know that in the maximum 
disorder case (T oo), the mobihty edge is \Ec\ ~ 3.6i. 
This energy is very close to the effective band edge and 
only less than 0.5% of the total states are localized. For 
lower temperatures we expect Ec to be closer to the effec- 
tive band edge and the number of localized states should 
be even smaller. We are interested in values of a: >0.1 for 
which the chemical potential is much higher than Ec and 
for these electron concentrations we expect the system 
to be metallic at any T. The DE model is an insulator 
only at very small values of the electron concentration, 



X <0.05 and high temperatures. With respect to the 
elastic mean free path, point b), we find that the con- 
ductivity of the system decreases continuously with T, 
until it reaches the T ^ oo limit (note that G ^ t). 
On the contrary if I is bigger than the cluster size the 
transport is ballistic, the conductance is determined by 
the cluster size, and increases as 7V^. When ^ is shorter 
than N the transport is diffusive and the conductance is 
proportional to N . 

In Fig. 3 we plot the conductance divided by N for 
x=0.5 and two values of the cluster size, iV=16 and 
iV=24. For low temperatures the core spins disorder is 
very weak and the elastic mean free path is larger than 
the cluster size. In this regime the transport is ballistic 
and the conductivity increases linearly with N . For large 
values of the temperature, G/N is independent of N . 
This is because the elastic mean free path is shorter than 
the cluster size and G/N is the electrical conductivity 
in units of the system lattice parameter. For T > 0.12t 
the values of G/N obtained using TV = 16 and iV = 24 
coincide, this implies that for these temperatures G/N 
is the conductivity. For x = 0.5 the magnetic critical 
temperature is Tc ~ 0.15^, and we can see in Fig. 3 that 
the conductivity is a smooth function of T, practically 
constant, near Tc- 

In Fig. 4 we present the same quantity G/N as a func- 
tion of T for X ~ 0.2. Also for this electron concentration 
the conductivity near Tc is a smooth function of T and 
the system is metallic. We have obtained similar results 
for different values of a; > 0.1. The values for Tc are 
taken from reference 1201. 
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FIG. 5. Conductivity, evaluated at Tc, as a function of x. 
The dashed line is only a guide to the eye. 
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In Fig. 5 we plot as a function of x, the value of the 
conductivity in lattice units evaluated at Tc- Using the 
value of AA for the lattice parameter, we obtain a resis- 
tivity p ^ O.OOlfi — cm at X = 0.2 near Tc. This resis- 
tivity is around ten times smaller than the experimental 
ones, which also have an insulator-like dependence with 
T. This is a clear indication that it is necessary to add 
other terms to the double exchange Hamiltonian, in order 
to explain the occurrence of the metal insulator transition 
at temperatures near the magnetic critical temperature 



V. SUMMARY 

We have calculated the temperature dependence of the 
electrical conductance of the double exchange Hamilto- 
nian. The Kubo formula has been use for the calculation 
of the conductance. Conductance is defined as the aver- 
age conductance over an statistical ensemble of clusters. 
These clusters are obtained by performing Monte Carlo 
simulations on the classical spin orientation of the dou- 
ble exchange Hamiltonian. The calculations have been 
done for different electron concentrations. We find that 
the system is metallic at all temperatures and, contrary 
to the mean field calculations, we do not observe any 
change in the temperature dependence of the resistivity 
near the magnetical critical temperature. Near Tc the 
resistivity we obtain is around ten times smaller than 
the experimental value. We conclude that the double ex- 
change model is not able to explain the metal to insulator 
transition which experimentally occurs at temperatures 
near the magnetic critical temperature. Other effects not 
included in the DE Hamiltonian, as electron-electron in- 
teraction or electron- phonon interaction, are needed in 
order to understand the electrical behavior of the oxide 
Mn. 
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